Three body recombination of ultracold dipoles to weakly bound dimers. 
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We use universality in two-body dipolar physics to study three-body recombination. We present 
results for the universal structure of weakly bound two-dipole states that depend only on the s-wave 
scattering length (a). We study threshold three-body recombination rates into weakly-bound dimer 
states as a function of the scattering length. A Fermi Golden rule analysis is used to estimate 
rates for different events mediated by the dipole-dipole interaction and a phenomenological contact 
interaction. The three-body recombination rate in the limit where a 3> -D contains terms which scale 
as a*, a^D^ and D'^, where D is the dipolar length. When a <^ D, the three-boby recombination 
rate scales as _D*. 
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Universality has become increasingly important and 
apllicable in ultracold physics. Two remarkable exam- 
ples are: strongly interacting fermions and few- 
body systems with large scattering lengths. Three-body 
recombination, A -\- A -\- A A2 -\- A -\- KE, is of spe- 
cial interest in an ultracold gas. It usually is the limiting 
process determining the lifetime of the gas in the absence 
of twp body losses. Three-body recombination has also 
proven to be an important probe in understanding the 
basic quantum mechanical nature of universal few-body 
systems Q. The experimental observation of universal 
3- and 4-body loss features 0-01 1 associated with a class 
of three body-bound state predicted by Efimov Q, is 
one of the major success stories in this field. These few 
body complexes have found great utility because of the 
control over the 2 body scattering length with a mag- 
netic Fano-Feshbach resonance @. In a magnetic field, 
the last bound state can be tuned through the scatter- 
ing threshold giving complete control over the scattering 
length. When a is much greater than the range of the 
short range interaction, many properties become inde- 
pendent of the details of the two-body interaction. For 
example, the three-body recombination takes on univer- 
sal scaling behavior of in this regime [loj . 

The powerful idea of universality can be exploited to 
study complex systems. Now ultracold polar molecules 
have emerged [ill, [13, these systems will have a 
new type of collisional control [l3| . This control emerges 
in the form of the dipole-dipole interaction, Vdd = 
— 3(z • f)'^]/r^ where r is the inter-particle separa- 
tion vector and d is the induced dipole moment along 
an external field axis (z). The collisional control can be 
achieved through the induced dipole moment d which is 
determined by the strength of the electric field. As the 
electric field is increased, there is a series of s-wave dom- 
inated resonances |l4| - [la |. These resonances occur when 
the system gains a weakly-bound long-range state. The 
short-range bound states structure of the interaction is 
largely independent of the external fields strength. These 



threshold resonances could be used to control scattering 
length, but the importance of these resonances will go 
further. The interaction is a long range and anisotropic 
which will mediate fundamentally novel quantum behav- 
ior. What is striking about this interaction is that the 
length scale describing it can be incredibly large, many 
orders of magnitude larger than the short range interac- 
tion. The dipolar length scale is 13 = fi2bd^/h'^ and ^26 is 
the reduced mass of a two body system. This large length 
scale can be used to create a highly correlated system. 

The two-dipole scattering cross section has been 
demonstrated to be universal [l^ [13, [ll] ■ In this work 
we extend the concept of universality in two-body dipo- 
lar physics and apply it to three-body recombination. We 
present results showing the universal structure of weakly 
bound two-dipole states that depend only on a. Using 
this, we then study threshold three-body recombination 
rates into weakly-bound dimer states as a function of a. 

To start, we study the universality of the dipolar 
threshold resonances. These universal dipolar collisions 
depend only on the s-wave scattering length. To tune the 
scattering length, we vary the short range position of a 
hard wall boundary condition [l^. As the inner wall's 
postition is decreased the system has a series of thresh- 
old resonances; in a real system this mimics increasing 
the electric field and the molecule's polarization. In this 
work we only study threshold resonances, the wide s-wave 
dominated resonances. There are other resonances which 
are not s-wave dominated, but these will be narrow and 
rare [Tsj . 

Examples of universal two-body scattering are well 
known, for instance when the scattering length is large 
and positive, large dimers form with binding energies 
Eb = h? I {2^2ba^)- Weakly bound two-body dipolar 
states near an s-wave resonance (a 3> D) also exhibit this 
universal property. The binding energy for these states 
is shown in Fig. [T] (a) in units of the dipolar energy 
Ed ~ h?/fJ,2bE>^ as a function of a/D. We have plotted 
the binding energies for three different resonances, the 
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first (black circles), second (blue diamonds) and eighth 
(red squares) resonance. As expected the three reso- 
nances agree very well in the large scattering length limit, 
but surprisingly the agreement is very good even when 
a<D. 

To further illustrate universal features, we look at 
the partial wave populations as a function of a. We 
have numerically obtained the molecular wave function: 
Tpd{ri2) = J2i Yio 12) fdi'^ 12) where I is the partial wave. 
We require that J2i^i = ^ where Ni = (/^|/^) is the par- 
tial wave population. Figure[T](b) shows the s- (solid), d- 
(dashed), and g-wave (dotted) populations for the first 
(black) and second (blue) resonances as a function of 
a/ D. This figure shows that the partial wave populations 
are universal. Additionally it shows that these molecules 
are mostly s- and d-wave. For large a, the s-wave con- 
tribution dominates and is near unity for large a and ipd 
takes on the universal form ripd oc exp [~r /a). As a is de- 
creased the d-wave contribution becomes significant and 
reaches ~ 40%. Even for very small a, the g-wave, i.e. 
I = 4, contribution is at most 2%. The dipolar interac- 
tion conserves parity and this selection rule prevents an 
s-wavc channel coupling to any odd partial wave. This 
implies that the present theory applies to only bosonic 
or distinguishable dipoles. 

There is of course a caveat to universality. The dipole 
length, D, must be much larger than any short range 
length scale, such as the depolarization, "dipole flip" 
length scale which is typically 100 ao where oq is the 
Bohr radius [l^. At this point the short range physics 
will become important and the structure of the bound 
state will be system dependent. However, these results 
are widely applicable because the significant size of the 
dipolar length scale. It can easily reach 10'^ ag, some can 
even exceed lO^ao for heavy, very polar molecules, such 
as LiCs. 

To find the three-body recombination cross section we 
begin by considering three distinguishable dipoles and 
use a simple Fermi's golden rule based approach which 
begins with the differential rate given by 
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''5{E,-Ef). (1) 



Here the initial state, ^i+i+i, is a three-body, box- 
normalized plane-wave in Jacobi coordinates. The fi- 
nal state, ^'2+1, consists of a weakly bound dimer and 
a box-normalized dipole-dimer plane wave. Because the 
dipole-dipole interaction has no diagonal s-wave contri- 
bution, we consider a modified two-body potential in Eq. 
([!]) which incorporates the zero-range Fermi-pseudo po- 
tential to account for s-wavc scattering, i.e. 
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where a is the positive s-wavc scattering length near a 
dipole-dipole threshold scattering resonance and 5 is the 
standard Dirac J-function. Summing over all final-state 
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FIG. 1: (Color Online) Universal properties of 2 body dipo- 
lar scattering, (a) The binding energy (symbols with curves) 
of the dipolar dimer as a function of the s-wave scattering 
length, for the first (black circles), second (blue diamonds), 
and eighth (red square) resonance, (b) The partial wave 
molecular populations are shown for the s-wave (soUd), d- 
wave (dashed) and g-wave (dotted) populations for the first 
(black) and second (blue) resonance. 



energies and dividing by the total incident particle fiux 
yields the total differential recombination cross section 
for three distinguishable dipoles [20j . 



d& 



dist 



d?p,d?p2rd{pi)e-'^''-'' (3) 



X y (fi3)e'(^i-'^i+^^-p^) 



Here pi is a mass scaled Jacobi vector given by 
(^■1 - r2) 
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In Eq. ([3]) ^2 is the outgoing dipole-dimer wave number 
with magnitude = \/2p,4^B (Emc + Eh), fci(2) is the in- 
coming wave number associated with the first (second) 
Jacobi vector 



k = ^/kl + k2 = \/2p.4BEinc is the mag- 
nitude of the initial three-body plane wave with total 
incident energy Ei^c- It should be noted that only a sin- 
gle two-body interaction is considered in Eq. This is 
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the simplest matrix element that connects the outgoing 
dipole-dimer state to the incoming plane wave state. 

Expanding the plane waves and dimer-wave function in 
Eq. ([3|) in spherical harmonics and integrating over the 
angular degrees of freedom yields the total distinguish- 
able particle cross section: 
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Here Jl is a spherical Bessel funtion, L and I are the 
initial angular momentum quantum numbers associated 
with the first and second Jacobi vectors respectively, L' 
is the angular momentum of the outgoing dipole-dimer 
system, and I' is the angular momentum of the dimcr. 
We note that we have moved to standard Jacobi co- 
ordinates in these integrals for ease of calculations, i.e. 
R = \/m/Mi2,3P2 and r = \/^jJ^2Pi- The ^Lo^io in Eq. 
(HI is due to the strong suppression all incoming angu- 
lar momenta greater than zero in s-wave scattering. The 
expectation value {L'V |Vd \ Ll)^ is only over aU an- 
gular coordinates. The cross section in Eq. 2] was found 
by summing over all final outgoing momenta, k'2 , and av- 
eraging over all incoming directions | fci , ^2 1 as would be 

appropriate for a gas phase experiment. 

The interaction in Eq. (f5|) is most easily evaluated with 
a biharmonic expansion |2l| . and it can be written as: 
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where r> (r<) is the greater (lesser) of ri and r2 and 
[Ykm Yk+2m']20 IS a rank 2 tensor formed by the two 
spherical harmonics. This can then be evaluated by con- 
sidering channels made up of two sets of partial wave, 
one for both fi and f2- 

With this long range interaction one might expect that 
all partial waves contribute equally to the recombination 
process as is the case in two-dipole elastic scattering. This 
does not occur for the three-body inelastic scattering pro- 
cess. As long as the incoming particles arc in the thresh- 
old regime, -Emc ^ Eb, all non-s-wave incoming channels 
are suppressed. This means that the final sum can be 
taken with zero initial angular momentum, L = I = 0, 
greatly simplifying the expression. 

The final identical boson cross section 0-*'''="* is found 
simply by multiplying by the Bosc enhancement factor, 
3!, and restricting the sum in Eq. ^ to outg oing chan- 
nels with the appropriate bosonic symmetry [2^. The 




FIG. 2: (Color Online) The three-body recombination rate 
of three identical bosonic dipoles is shown in dipole units as 
a function of the scattering length. The total recombination 
rate (solid curves) is universal and independent of the reso- 
nance used. The sS (dashed red), dS (dotted red) and higher 
order (dot-dashed red curves) partial wave contributions are 
shown as well. 



incoming channel is restricted to s-wave only and al- 
ready obeys the appropriate symmetry. The recombina- 
tion rate can be extracted from the inelastic cross section 
by multiplying by the characteristic incoming velocity 
V = hk/^iiB [lfl[2i: 
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In general the integrals in Eq. ([S]) must be per- 
formed using numerically tabulated dimer wave func- 
tions (r), but in the case where a ^ D Fig. [T] (b) 
shows that the dipolar wave function is dominated by 
the s-wave component. Furthermore the wave function 
has the universal form for a weakly bound s-wave dimer: 
V'd (r) e^''/''/(rV27ra). Inserting this into Eq. (O 
yields the large a scaling of the recombination rate: 
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The leading order term in this is exactly that found 
in Ref. j22l [23| for identical bosons with short-range s- 
wave interactions with the appropriate scattering matrix 
element set to (ka) . The second term is due entirely to 
the presence of the dipolar interaction and corresponds to 
an outgoing d-wave differential cross section. It should be 
noted that higher partial wave contributions in the dimer 
wave function are present at the a^D^ order as well, but 
this contribution is two orders of magnitude smaller. 

The results of the Fermi gold rule calculations are 
shown in Fig. [5] as a function of scattering length for 
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the first resonance in Fig. [T] The total rates have the ex- 
pected scaling at large a due to the contact interaction, 
but as a/D is decreased they eventually flatten out and 
when D > a the rate is entirely controlled by the dipolar 
interaction, scahng as D"^: Kf""-^ 0.079hD'^/ fj.2b- We 
have compared the recombination rates for several two- 
body resonances and have found nearly perfect agree- 
ment, to within ~ 1%, indicating that the three-body 
recombination behavior found here is universal. 

The dipolar interaction directly couples the incident 
channel to final states with L' = 0,2,4, ... and a molec- 
ular state with I' = 0,2,.... The largest term in the 
a > D regime, I'L' = sS (dashed red curve), comes 
from the contact interaction and scales as a^. The next 
largest contribution from the sD channel (dotted red 
curve) scales as D^a?. Additional contributions from 
higher order I'L' channels arc shown (dot-dashed red), 
but contribute little in the large a regime. However in 
the a < D regime, the dominant contribution come from 
the I'L' = dS and dG terms, the first of which is only 
weakly dependent on the scattering length. As one might 
expect there are two universal regimes that appear. The 
first regime appears where the long-range dipole inter- 
actions is dominant and the scattering length is small, 
i.e. a <C £>, where is independent of a and scales as 
_D^. The second regime appears when a is the dominant 
length scale, a ^ D, and the s-wave contribution is 
dominant over all others. 

The Fermi's golden rule based approach is only capable 
of describing the scaling behavior of three-body recom- 
bination, and does not incorporate more complex three- 
body correlations. The recombination rate can have a 
variety of interesting resonant behaviors on top of the 
envelope behavior presented here. For instance we ex- 



pect that in the a ^ D regime, a geometrically spaced 
set of Efimov minima are likely to appear in the I'L' = sS 
outgoing channel [3l ■ 

We have shown the universal structure of weakly 
bound two-dipole states that depend only on the scat- 
tering length of the system, shown in Fig. [T] This has 
allowed us to obtain threshold three-body recombination 
rates into weakly-bound dimer states as a function of a. 
We used a Fermi Golden rule analysis to estimate the 
contributions to the rate from different events mediated 
by the dipolc-dipole interaction and a contact interac- 
tion. These rates are shown in Fig. [2] The individ- 
ual terms have many different scaling behaviors. When 
a > D, we find the dominant term scales as a"*. When 
D > a the recombination rate scales as D'^. The case of 
negative scattering length is not addressed in this work, 
but it is not unreasonable to expect that the envelope 
behavior will be similar to that presented here. More 
detailed three-body calculations which include extra res- 
onance structures are the subject of future work. The 
|a| ^ D regime might also hold a variety of interesting 
phenomena such as universal thrce-dipole Efimov states. 
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